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Abstract 

The cosmological evolution of neutrino energy densities driven by cosmon-type field equations 
is introduced assuming that active and sterile neutrinos are intrinsically connected by cosmon 
fields through the seesaw mechanism. Interpreting sterile neutrinos as dark matter adiabatically 
coupled with dark energy results in a natural decoupling of (active) mass varying neutrino (MaVaN) 
equations. Identifying the dimensionless scale of the seesaw mechanism, m/M, with a power of 
the cosmological scale factor, a, allows for embedding the resulting masses into the generalized 
Chaplygin gas (GCG) scenario for the dark sector. Without additional assumptions, our findings 
establish a precise connection among three distinct frameworks: the cosmon field dynamics for 
MaVaN's, the seesaw mechanism for dynamical mass generation and the GCG scenario. Our results 
also corroborate with previous assertions that mass varying particles can be the right responsible 
for the stability issue and for the cosmic acceleration of the universe. 
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The mass varying mechanism has often been suggested as a cosmological clock for the 
time evolution of the dynamical dark energy [l, 2]. Theories that describe growing matter 
in such a context can possesses an adiabatic regime in which the quintessence field always 
sits at the minimum of its effective potential, which is set by the local mass varying dark 
matter density. That is the prerogative of the so-called stationary condition of mass varying 
mechanisms [3H5|. As a direct consequence of this exotic interaction, the relic particle mass 
is generated from the vacuum expectation value of the scalar field and becomes linked to 
its dynamics. This coupled fluid is usually identified either as dark energy plus neutrinos, 
or dark energy plus dark matter 12]. In mass varying neutrino (MaVaN) scenarios, the 
time evolution of a cosmological scalar field can be interrupted by an increasing mass of 
the neutrinos. Some approaches predict a transition from a cosmological scaling solution 
with dynamical dark energy at early times to a cosmological constant dominated universe 



at late times 



a a. 



In a particle physics realization of such a growing matter scenario, 



MaVaN's act as a trigger for the scalar field dynamics, giving an explanation to the "why 



now problem" [2j,l2j. The interface between MaVaN's and mass varying dark matter can also 
provide interesting relations among present values of neutrino masses, dark energy density, 



a. 



their corresponding equations of state and the cosmological stability conditions 

In this letter we provide some convincing constraints among MaVaN's, dynamical dark 
matter, and dark energy, through a single degree of freedom: a dynamical scalar field, 4>(a), 
or more explicitly, the cosmological scale factor, a. Through a natural connection with the 
generalized Chaplygin gas (GCG) js, 9|, we show that the intermediacy of the dark sector in 
teractions done by 0(a) can be exclusively relegated the scale factor, a, and correlated scales 
cosmologically driven by the GCG equation of state. The first prerogative of our approach is 
that the cosmological evolution of active and sterile neutrino energy densities are driven by 
symmetrically coupled cosmon-type field equations [ll|. Active and sterile neutrino states 
are thus intrinsically connected by cosmon fields through the seesaw mechanism for mass 
generation 10(. Active neutrinos presumably have tiny masses when compared to charged 
fermions, or to supermassive sterile neutrinos. Such a smallness is understood in terms of the 
symmetry of the standard model (SM) of electroweak (EW) interactions amended by seesaw 
mechanisms. As suggested by previous issues jl|, the sterile neutrino mass, M. = Ai(M(4>)), 
exhibits dynamical behaviour driven by the dark energy cosmon field, <fi. Interpreting sterile 
neutrinos as the aforementioned dark matter adiabatically coupled with dark energy results 
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in active neutrinos with dynamical mass, i. e. MaVaN's. They turn consistently into 
a perturbative effect driven by the abovementioned cosmon fields. It results in a precise 
description of the dark sector energy components: dark energy dark matter (as sterile 
neutrinos) and (active) neutrinos. We consider that cosmon fields interconnect the dark 
sector and MaVaN components through the seesaw mechanism. 

Dark matter is most often not considered in the formulation of the MaVaN models; how- 
ever, the possibility of treating dark energy and dark matter in a unified scheme naturally 
offers this possibility. Therefore, through the second prerogative of our approach, we assume 
that the mass varying behaviour, the cosmological dependence of the energy components, 
and the conditions for cosmological stability are all guaranteed by a quite interesting con- 
nection between the mechanism for mass generation and the GCG. Its dependence with 
the scale factor indicates that the GCG can be interpreted as an entangled admixture of 
dark matter and dark energy. We shall verify that identifying the dimensionless scale of the 
seesaw mechanism, m/M, with a power of the cosmological scale factor, a, allows for embed- 
ding the resulting dynamical masses into the GCG scenario for the dark sector. Thinking 
on an effective model approach, it eliminates the intermediacy of dark sector interactions 
by a scalar field, in favor of simple scale dependencies. Without additional constraints, our 
findings establish a precise connection among three distinct frameworks: the cosmon field 
dynamics, the seesaw mechanism for dynamical mass generation and the GCG equation of 
state. 

The Lagrangian densities of active (A) and sterile (S) components in terms of matter 
fields, ip A ,s, can be written as 

£a,S = i4>A,sl^*pA,S + k A , S 4> A , S tp A , S , (1) 

where we have introduced two mass scales, k A = \i and k s = Ai, which respectively designate 
the tiny and large masses predicted by the seesaw mechanism, 

-H = (M/2)[l - y/l + 4(m/M) 2 ] and M = (M/2)[l + y/l + 4(m/M) 2 ]. (2) 

m 
m M 

orthogonal basis of chiral left- and right-handed neutrinos, z/^ R , related with the matter 
fields, Va,s, by 

i> A = Vs 2 + l(sz/£ - u R ) and = V s 2 + l(z/£ + sv R ), (3) 



They correspond to the exact eigenvalues of the mass matrix 



written in the 



where we have introduced the the dimensionless scale s = \J Jl/M. It results in C A + C s = 
£-h + + £-\j R i where £^ give the kinetic contribution and C-^ lR gives the mass mixing 
terms 10|. The equivalence between the stress-energy tensor in chiral basis and matter 
field basis results in p A + p s = P£ + p R and p A + p s = p^ + p R . After introducing an 
auxiliary mass scale £ = \fp~JA = m, we provide the ingredients to define two energy scales, 
Pi — (Ps + Pa)/2 and p 2 = (p s — p A )/2, which cosmologically evolve driven by symmetrically 
coupled cosmon-type equations. In a FRW universe it corresponds to 

p 1)2 + 3H( Pl , 2 + p 1>2 ) - + = o (4) 

where if = a/a is the expansion rate of the universe, and the indices 1, 2 denote the 
symmetry under the interchange of p 1 by p 2 . The third term is the mass varying term and 
the last one is a mutual exchange term due to the coupling between states 1 and 2. p A and 
p s are evaluated like normal modes of p x and p 2 which are driven by symmetrically coupled 
cosmon-type equations. p 1 component can be identified in terms of matter fields while p 2 is 
just an auxiliary energy scale that evolves with cosmon dynamics. 
After combining Eqs. ((2]), and observing that, 

p A (a,£,s) = p A (a,£s) (5) 

and 

p s (a,^s) = p s {a,£/s) (6) 

it is easy to identify the effective dynamical masses explicitly related to the cosmological 
evolution of active and sterile energy densities, p A and p s , by means of decoupled equations, 

Pa + 3H{Pa+Pa) .^ = , (7) 



and 



A + Stf(* + *)-^^ = 0. (8) 



dm£ .„ /dlnu d\nM 

- '1/2) + 



from where one notices that 



d<p \ d(j> d<p 

and 

dins ( dlnp dlnA^ 



d0 V d(p d(f> 



The coupling between relic particles and the scalar field as described by Eqs. ©-(jSJ) are 
effective just for non-relativistic (NR) fluids. Since the strength of the coupling is suppressed 
by the relativistic increase of pressure, as long as particles become ultra- relativistic (UR), 
with T(a) = T /a ^> m(<f>(a\), the dark matter fluid, p s , and the neutrino fluid, p A , decouple 
and evolve adiabatically That is the case for ultra-relativistic (UR) regimes where 
~ |j oc (p — 3p) ~ 0. The mass varying mechanism expressed above translates the 
dependence of masses on the scalar field into a dynamical behaviour, i.e. p(a) and Ai(a). 

In some previous issues, we have already suggested that one could treat MaVaN's as a 
oerturbative component derived from some previously unperturbed adiabatic solution p^ 



12j. The above results provide a precise explanation for this. As one can observe, all 
the information from the dark sector (dark energy plus dark matter) acting on the (active) 
neutrino sector is carried out by the explicit dependence of p(4>) on <p. From the cosmological 
point of view, it results in matter fields (active and sterile) that evolve separately. It is no 
way the case for the coupled chiral eigenstates and if) R . At primordial times, when s 2 ~ 1, 
such mass eigenstates should be indistinguishable, and the chiral eigenstates symmetrically 
defined. At late times they turn into dark matter and (active) neutrinos, maintaining the 
correspondence among the flavour sectors. To proceed, one can notice that the cosmon 
framework provides the right connection between dark matter (sterile neutrino) and dark 
energy through the well-known cosmon field equation, 

p*+mp*+P*)+^^=o- (9) 



'nsl-i 



Eq. ([?]) for i — A and Eq. (JUJ) result in the adiabatic equation for the dark sector (p DS ,Pi 

p DS + 3H(p DS +p DS ) = 0, (10) 

with p DS = p^ + p s = p^ + p ± + p 2 and H 2 = p DS (with H in units of H and p DS in units of 
Pent = 3H 2 /8t[G). In spite of an intrinsic dependence on 0, the equation of motion for the 
dark sector is not modified by p A . The cosmological evolution of p A can be perturbatively 
computed through Eq. (JTj). 

The phenomenological consistency of this scenario can be quantified when one approx- 
imates p A and p s by the energy densities of a degenerate fermion gas (DFG) at different 
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relativistic regimes, 

p A (a) = (87T 2 )-Xa) 4 Ua)(2r/(a) 2 + l)^r]{a) 2 + 1 - arc sinh (77(a)) 



p s {a) = (87r 2 )- 1 /i(a) 4 s- 8 7^77(0) (27 Vr/(a) 2 + 1) vV^a) 2 + 1 - arc sinh (-fs 2 r](a)) (11) 

where 77(a) = q A /(ap(a)), g F is the Fermi momentum, and the relation between the fluid 
thermodynamic regimes is parameterized by the coefficient 7 = qs/q A - in particular, the 
DFG prescription is adequate for reproducing a consistent analytical transition between UR 
and NR regimes, and the inherent effects due to coupling dark matter (Ai (</>)) and dark 
energy (0) can be verified through Eqs. (JTJ)-(J9]) . 

In the Fig. [T]one observes the exact correspondence between the abovementioned energy 
density components and the "modified" scale parameter 7s 2 . For the NR limit of a DFG, 
we have p A /ps ~ 7" 3 s 2 . 

The characteristic size of the active neutrino masses involves appropriate combinations of 
dimensionless (Yukawa) couplings Y jy i. e. /x, ~ m 2 /M. with ~ Y 3 v (v ~ 2x 10 11 eV), and 
the consistency with the observed oscillations requires for the mass of at least one neutrino 
A*j ~ 0-05 eV. For Yj of the order one this implies an upper bound for the large mass 
scale M. < 10 23 eV, which results in s 2 > 10" 24 . Once the present value of the rate p A j ' p s is 
phenomenologically assumed as 0(1CT 2 ), a DFG of (active) neutrinos at least approximately 
in the NR regime leads to 7" 3 ~ 10 22 , i. e. two completely different momentum scales for 
active and sterile neutrinos, q s /<lA ~ 10~ 7 . Sterile neutrinos behave like ultra cold dark 
matter. The phenomenological value of p A /Ps ~ 0(1CT 2 ) is consistent with the ACDM 
cosmological model with (massive) cosmological neutrinos which transit from the UR to 



the NR regime in the CDM domination era js]. The super massive sterile component is 
supposed to transit to the NR regime at energies much larger than 1 MeV, i. e. earlier 
than the period where the relativistic particles: photons and electron/positrons were in 
equilibrium with neutrinos. Thus, ultra cold dark matter could have become NR much 
earlier than neutrino decoupling from the SM cosmological plasma. A basic understanding 
of the interaction rates of active neutrinos enables us to argue that active neutrinos were 
once kept in equilibrium with the rest of the cosmic plasma. Ultra cold dark matter lost 
contact with the cosmic plasma slightly before reach the NR regime, so it did not inherit 
any remaining cosmic plasma associated energy. The active neutrinos, then which did, are 
therefore hotter than the sterile ones. They could have been "re- heated" by SM interactions 



which remain active after CDM decoupling as well as photons decoupled from neutrinos 
were re-heated due to electron-proton recombination. It could result in the sterile/active 
momentum discrepancy given by q s llA ~ 10" 7 . Once the partial derivatives of energy 
densities with respect to masses in Eqs.([7HS]) can be written as 

dp A pA~3p A dp s p s -3p s 

op p DM M 

they establish when the coupling to the scalar field is effective. In particular, the condition 
given by Eq. fTTZj) is guaranteed by the DFG approach. 

In addition, if one assumes that dependence of the energy density on the scalar field 
could be parameterized by a \<j) 4 theory, or some type of quintessence potential the adiabatic 
evolution of the scalar field allows for obtaining the mass of the scalar field. Such predictions 
for are consistent with the lower bound on the mass derived from assumptions that the 
scalar field has evolved adiabatically since nucleosynthesis. The scalar field mass should 
then be greater than the Hubble constant at nucleosynthesis given by ~ 1CT 17 eV. However, 
the explanations for the adiabatic regime where the light scalar field have settled down to 
the minimum of its potential prior to nucleosynthesis are quite model dependent. It requires 
a huge fine-tuning to generate the different scalar mass scales and it should be embedded 
in a more encompassing model, such as the Minimal Supersymmetric SM with the addition 
of one singlet chiral superfield [13]. For a unified picture of dark energy and dark matter, 
the Higgs boson can be coupled to an additional singlet scalar field which we identify with 
a quintessence field [l4j. These predictions are consistent with the widely believed assertion 
that the SM of particle physics is actually only a low-energy approximation of a more 
fundamental underlying theory. Such an interplay with cosmology represents an important 
guideline to obtain insights on the nature of a more fundamental theory. It concerns, for 
instance, a natural connection with the cosmological scenario of the GCG. The GCG model 
is characterized by an exotic equation of state |S|, [9] given by 

P = -A Po (^J, (13) 

which can be obtained from a generalized Born-Infeld action [8]. In any case, irrespective of 
its origin, several studies yield convincing evidence that the GCG scenario is a phenomeno- 
logically consistent approach to explain the accelerated expansion of the universe. Notice 
that for A s = 0, the GCG behaves always as matter whereas for A 3 = 1, it behaves always 



as a cosmological constant. This property makes the GCG model an interesting candidate 
for the unification of dark matter and dark energy, i. e. for the dark sector energy density, 
p DS , of our model. 

Inserting the above equation of state into the unperturbed energy conservation Eq. (|10p . 
one obtains through a straightforward integration |8j, 



Pds = Po 



A a + 



.4. 



l/(l+«) 



(14) 



3(i+a) 

One of the most striking features of the above equations is that the GCG energy density 
interpolates between a dust dominated phase, p DS oc a -3 , in the past, and a de-Sitter phase, 
p DS = —Pds, at late times. This evolution is controlled by the positive parameters, a and A s 
with < a < 1 . Of course, a = corresponds to the ACDM model and we are assuming 
that the GCG model has an underlying real scalar field js]. All the phenomenological issues 
of the GCG have been recently addressed in Ref . [l^| . 

Assuming the canonical parametrization of p DS and p DS in terms of a scalar field <p, 

1 



Pds = ^<f 2 + an d Vl 



2^ 



(15) 



and the Friedmann equation H 2 = p^ in natural units with Planck mass Mpi = 1 (i. e. with 
H in units of H and p D s in units of p Crit = 3H 2 /8nG), one obtains the effective dependence 
of the scalar field (p on a given by 

Vl + a^AJ(l - A.) - 



<p(a) 



2/3 



(16) 



Y 1 + a^A s /(l - A s ) + 1_ 

where (3 = 3(a+l)/2. From this point, the dynamics of the GCG scalar field is parameterized 
by Eq. f|T6l) . Thus we can consider the simplest version of the seesaw mechanism for which 
the Yukawa coupling between a light scalar field and a single neutrino flavour prescribes a 
linear dependence of the mass scale m on if, i. e. m((p) ~ ip. By observing that the logarithm 
of the squared scale parameter s 2 has an analytical structure similar to that of <p(a), 

1 



ln(s 2 ) = In (p/M) = In 



• v /l + 4(m/M) 2 



(17) 



l + 4(m/M) 2 + l 

we can rewrite the auxiliary scale m/M in terms of the GCG parameters A s and (3, and of 



the scale factor a, as 



m 



m(a) 

M{a) ~ M 



A. 



I — A, 



ina 



(18) 
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with the arbitrary constant k obviously introduced for phenomenological reasons (k disap- 
pears if one redefines A a and p in Eq. (T14]) ). After simple mathematical manipulations, we 
obtain 

In (s 2 ) = In (n/M) = -2flip( K a), (19) 

Adjusting the Yukawa coupling to the GCC scalar field parameters leads to m(<j)(a)) = 
Y<j)(a) ~ ip(na). Eqs. ffTTj) -( fT9"j) result in the following simplified forms for active and sterile 
neutrino masses, 

= ip exp (— flip) and M.{ip) = ipexp (+fl<p). (20) 

where <p = ip(na). Since the scalar field constrained by the GCG dynamics decreases with 
a, for fl > 0, active neutrino masses, /z, increase with a, and sterile neutrino masses, Ai, 
decreases with a. Such a behaviour sets an obvious analytical dependence of the scale pa- 
rameter, s, on ip given by an exponential function: s = exp (— flip (a)). Such an exponential 
divergency between two mass scales is consistent with the naturalness of exponential po- 
tentials for several classes of quintessence models J^, [lsj. In addition, MaVaN's can be 
perturbatively implemented through Eq. <^j. As in some previous issues, growing neutrinos 
essentially stops the cosmological evolution of the scalar field and triggers an accelerated 
expansion of the universe. It is evident that our findings are fairly general, as well as 
independent of the choice of the equation of state, and the form of M(4>) and modifica- 
tions on the equation of state for dark energy can lead to quite different scenarios. For 
instance, it also admits static dark matter, which should result in neutrino masses given by 
/i(a) oc exp (-flp(Ka)) = exp (— flY 0(a)). Turning back to the possibility of adiabatic insta- 
bilities, they were previously pointed out in cosmological scenarios in a context of MaVaN 
models for dark energy. The essential ingredient of this class of models is that neutrino mass 
depends on the dynamical dark energy and grows in the course of the cosmological evolu- 
tion. The tiny neutrino mass and the recent accelerative era are twinned together through 
a scalar field coupling. In the adiabatic regime, these models face catastrophic instabilities 
on small scales characterized by a negative squared speed of sound of the effective coupled 
fluid. Such instabilities would give rise to exponential growth of small perturbations over 
the background uniform fluid [6|. The natural interpretation of this is that the Universe 
becomes inhomogeneous with neutrino overdensities subject to nonlinear fluctuations which 
eventually collapses into compact localized regions. 

9 



The mechanism that we have prescribed naturally leads to a fast increase of the neutrino 
masses, which results in a quite model dependent vanishing squared speed of sound at the 
present epoch. The squared speed of sound is kept positive during the entire cosmological 
evolution The GCG, as the cosmological background fluid, is consistent with all these results. 
In particular, the dynamical mass is due to the time evolution of a cosmon field that coincides 
with the dynamics prescribed by a GCG universe. These scenario is by no means the 
only possibility. However, without any additional assumptions, our findings provide an 
attractive and convincing explanation for the why now problem through a confluence among 
three independent frameworks: the cosmon field dynamics, the seesaw mechanism for mass 
generation and the cosmological GCG scenario. For active neutrinos, an increase of the 
mass n by a factor 10 6 approximately corresponds to a decrease of the sterile neutrino mass 
M. from the Planck scale to 10 13 GeV. It would largely be sufficient to load to dramatic 
consequences for cosmology, even in the present cosmological epoch. As soon as neutrinos 
become non-relativistic, their coupling to the cosmon triggers an effective stop (or substantial 
slowing) of the evolution of the cosmon. Naturally our results deserve a deeper analysis in 
which concerns the fine-tuning of the phenomenological masses involved. 

The main features of our model can be summarized in the following terms, i) The cosmo- 
logical evolution of neutrino energy densities is driven by symmetrically coupled cosmon-type 
field equations that result in unmixed equations for active and sterile neutrino states intrin- 
sically connected through the seesaw mechanism for mass generation, ii) The possibility of 
treating dark energy and dark matter in the GCG unified scheme naturally couples dark 
matter with MaVaN's. Interpreting sterile neutrinos as dark matter adiabatically coupled 
with dark energy provides the sufficient conditions to implement such unified picture in Ma- 
VaN formulation. We have found the constraints imposed by the seesaw mechanism in order 
to establish the unique analytical connection to the canonical formulation of the GCG, in 
terms of a real scalar field, iii) The dynamics of the mutual coupling among neutrinos, dark 
matter and dark energy are driven by only one degree of freedom: the scalar field, (pa. As is 
well known, the Higgs sector and the neutrino sector are possibly the only ones where one 
can couple a new standard model (SM) singlet without upsetting the known phenomenol- 
ogy. Unfortunately, we know to much about the Higgs of the SM seesaw mechanism, i. e. 
the constraint on its interactions and properties. It drops any pretension of classifying the 
generic scalar field of our model. Since the GCG is formulated in terms of its equation 
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of state that results in an explicit dependence on the scale factor, a, the introduction of 
the scales s and £ allows for eliminating the intermediacy of the scalar field, in favor of 
assuming an explicit dependence on a, that is p, p(4>(a)) as p, p(a). Due to the connection 
that we established between the GCG and the seesaw masses, our approach can also result 
in an effective model for MaVaN's coupling with the dark sector. Eliminating the coupling 
action of the scalar field, in this case, the dynamics is driven by dimensionless scales a, s(a), 
and £(a)/M(a). To conclude, our unifying picture predicts the same features of standard 
quintessence models in mass varying matter scenarios which, however, falls through more 
ambitious unifying schemes. Our results reproduce the quintessence exponential laws be- 
sides corroborating with previous assertions that mass varying particles can be the right 
responsible for the stability issue and for the cosmic acceleration of the universe. Acknowl- 
edgments: This work was supported by the Brazilian Agencies FAPESP (grant 08/50671-0) 
and CNPq (grant 300627/2007-6). 
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FIG. 1: The exact correspondence between the energy density components and the "modified" 
seesaw scale parameter 7s 2 . Sterile and active neutrinos are assumed to behave like a DFG. For 
the NR limit, we have Pa/ Ps ~ 7~ 3 s 2 . 
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FIG. 2: Growing neutrino mass Qu(a) / /j, ) , decreasing dark matter mass (Ai(a) / M Q ) and neutrino 
energy density (p u (a) / p u0 ) in dependence on scale factor a. They are respectively normalized by 
unitary quantities of fi , M and p u0 at present (a = 1). The results are obtained for different 
GCG scenario with A s = 4/5 and parameters a = 1 (solid line), 1/2 (dot line), 1/4 (dash dot line), 
and 1/8 (dash line). 
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